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Summary 

An algorithm suitable for a mini -computer was 
developed for finding the natural frequencies and 
mode shapes of a planetary gear system which has 
unequal stiffnesses between the sun/planet and 
planet/ring gear meshes. Node shapes are re- 
presented In the form of graphical computer output 
that Illustrates the lateral and rotational motion 
of the three co-axial gears and the planet gears. 
This procedure permits the analysis of gear trains 
utilizing non-uniform mesh conditions and user 
specified masses, stiffnesses, and boundary condi- 
tions* Numerical integration of the equations of 
motion for a three planet system was performed. 
Excellent agreement was obtained between natural 
frequencies calculated assuming constant mesh 
stiffnesses and resonant frequencies of the dynamic 
system. The algorithm can be used for determining 
the time step and total Integration times In the 
numerical solution of the differential equations of 
motion. The methodology is a quick and efHclent 
means of Identifying operating resonance regions 
and analytically tuning the drive train by examining 
the effect of parameters such as planet and sun gear 
support stiffnesses, input/output shaft stiffnesses, 
and gear mesh stiffnesses. The algorithm is adapt- 
able for use In a Computer-Aided Design (CAD) 
environment. 

Nomenclature 

Fx force in x-direction. N (lb) 

Fy force in y-direction, N (lb) 

planet index 

K translational support stiffness, N/m (Ib/in) 
KPR planet/ring mesh stiffness, N/m (Ib/in) 

KSP sun/planet mesh stiffnesses, N/m (Ib/in) 

KT torsional support stiffness, N/m (Ib/in) 

R base circle radius, m (in) 

T torque, N m (lb in) 

o angle of line of action for sun/planet 
mesh (see Fig. 4) 

0 planet reference angle (see Fig. 4) 

Y angle of line of action for planet/ring 

mesh (see Fig. 6) 

a gear pressure angle (see Fig. 4) 


August is presently employed by Life 
Systems, Inc., Cleveland, Ohio. 


Introduction 

Many advanced technology applications In power 
transmissions, such as aircraft applications, have 
a high power to weight ratio requirement. Planetary 
gear trains, rather than parallel shaft transmis- 
sions. are widely used for these applications 
because of a considerable savings In weight, space 
and drive train efficiency. These advantages stem 
from the fact that multiple planets allow the load 
to be transmitted by several tooth contacts. For a 
planetary arrangement with three planets, each 
tooth engagement of the sun gear must carry only 
one third of the total static load. Consequently, 
the dimension of the sun gear may be considerably 
reduced from that of the pinion of a parallel shaft 
gear designed to transmit an equivalent torque. 
Furthermore, the co-axial arrangement of driven and 
driving shafts allows a compact layout. 

The analytical determination of the natural 
frequencies of planetary gear trains Is an Important 
design criteria for three reasons. First, In addi- 
tion to considering the input/output speeds of the 
transmission, the power transmitted and tooth loads 
and stresses, the designer must take care to assure 
that he does not excite a natural frequency within 
the power transmission system. This may lead to a 
fatigue problem and eventual premature failure of a 
drive train component. 

Secondly this algorithm provides the designer 
with a tool with which he can analytically **tune'* 
his drive train. He may interactively modify impor- 
tant drive train parameters such as the stiffnesses 
of the planet-ring mesh, planet support. Input/ 
output shaft and transmission housing. 

A third reason, and the one that Initiated the 
research documented here, is to provide a basis for 
selecting a time step for the numerical Integration 
of the equations of motion of the planetary gear 
train. This integration is important in determining 
the dynamic loads, and stresses. Judicious selec- 
tion of the time step assures the stability of the 
integration process. 

Several investigators have considered the 
natural frequencies of planetary gear trains having 
equally loaded/epual stiffness planet meshes. 
Cunllffe, et al.^ was successful In correlating 
results from an experimental test rig and analog 
computer analyses to analytical determinations. 

Using these approaches, they modeled simple plane- 
tary gear trains as well as some complex designs 
considering linear and nonlinear mesh stiffness. 
Their Investigation was limited to a star gear 
system with a stationary carrier and a rotating 
annulus. This resulted In a system having 13 kine- 
matic degrees of freedom. 

Botman^ extended this work to an eighteen 
degree of freedom model that included the effects 
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of planet pin stiffness and rotation of the carrier. 
It 1$ the objectives of the analysis reported 
herein to (a) extend Botman's ivork to Include un>* 
equal aiesh stiffness between the three Planets In 
the sun*planet and/or planet-rlnq meshes* (b) to 
analyze the effect of a thin tooth In one of the 
components and (c) to simulate unequal stiffnesses 
due to alternation between one and two pairs of 
teeth In contact. 

Dynamic Nodel 

Background 

This analysis Is 11m ted to the free In-plane 
vibration (undamped) of a single-stage, multiple- 
planet system such as the one illustrated In Fig. 1. 
N’hlle all of the figures Illustrate a three-planet 
system, the method Is directly applicable to sys- 
tems with any number of planets. Components are 
treated as lumped-parameter rigid bodies. The gen- 
eralized Jacobi eigenvalue routine used requires 
symmetrical Inertia and stiffness matrices, there- 
fore, a global coordinate system Is requlr^. The 
general model shown In Figs. 2 and 3 Involves 
18 degrees of freedom: three for each component 
(one rotational, two translational). The components 
Include sun gear, ring gear (annulus), carrier, and 
three planets. Simpler models Involving fewer 
degrees of freedom can be simulated by suppressing 
unwanted modes with extremely high values of the 
support stiffness parameter. 

Each member 1$ modeled as supported by a tor- 
sional spring of stiffness KT and a bearing sup- 
port of two equal springs of stiffness K located 
along the X-Y coordinates axes. (Note: For the 
planets, KT Is zero.) Gear mesh Interactions 
(which are treated as parameters) are modeled as 
linear springs of constant stiffness acting tan- 
gentially at the respective base circles. The sup- 
port springs and mesh Interactions are shown In 
Fig. 3 (only one planet Is shown). Constant tooth 
stiffnesses permit a linear stiffness matrix and a 
manageable elgenroutine. Botman^ states that 
this approach predicts realistic frequencies and 
mode shapes. Calculation of dynamic tooth loads 
(not considered In this model) requires time vary- 
ing stiffnesses to account for the alternation 
between one and two pairs of teeth in contact at 
each mesh Interaction. 

Equations of Notion 

The equation of motion for an n-degree of 
freedom system in free vibration is: 

[M](q> ♦ [C]{q> ♦ [K](q> . {0> (1) 

where M is the Inertia (mass) matrix, C Is the 
damping matrix (assumed zero), K Is the stiffness 
matrix and q Is the vector of the generalized 
coordinates (a, x, y) for each component. For the 
lumped mass system here, the Inertia terms form a 
diagonal (uncoupled) matrix and the stiffness matrix 
Is symmetric. The stiffness terms shown In Table 1 
are developed by calculating the reaction (torque or 
force) resulting from a unit relative displacement, 
q. In each coordinate direction (a, x, y) for each 
component. Component Interactions are shown In 
Figs. 4 to 6. 


Assembly of Inertia and Stiffness Xatrlces 

The d1aqr>al Inertia matrix Is composed of the 
mass moment of Inertia term and two equal mass 
terms (representing the three degrees of freedom) 
for each component. The negative of each reaction 
term Is assembled Into a stiffness matrix by taking 
the terms In systematic order consistent with the 
or<ter used In the Inertia matrix. While any con- 
sistent ordering of terms may be used, a logical 
ordering with directly Interacting components 
placed adjacent to each other minimizes the band 
width of the elgenproblem. This produces faster 
convergence and Improved accuracy. The components 
are ordered: Sun, planet 1, planet 2, ..*, planet 
n, carrier, and ring. The three degrees of freedom 
for each component are taken In order: a, x, y. 
Thus, the stiffness matrix includes: 
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Each diagonal term of the stiffness matrix 
also includes the appropriate bearing support reac- 
tion. Since the stiffness matrix is symmetric, only 
the lower (or upper) triangle need be actually com- 
puted, the rest being loaded from symmetry. 

The natural frequencies and eigenvectors are 
calculated by the Jacobi subroutine [3] using a 
16-b1t-word mini -computer. Core requirements 
(including graphics software) fit within a 32 K 
word partition. 

Presentation of Mode Shapes 

Upon exit of the elgenroutine, the eigen- 
vectors are stored columnwise In the solution array. 
These represent normalized relative displacements 
of the components and are arranged In the same 
order as Is used In assembly of the Inertia and 
stiffness matrices: Motion of sun, planets no. 1 
through n, carrier, and ring. For each component, 
displacements are given In the order rotation, 
X-translatlon, and Y-translatlon. 

To aid Interpretation of planet motion, local 
coordinate motion (radial and tangential) may also 
be computed for the planets. (Note: The Global 
coordinates have been chosen such that planet no. 1 
is at 0 and. therefore. Its radial motion Is along 
the X-axis and Its tangential motion Is along the 
Y-axis: Local and global coordinates are identical 
for planet no. 1.) 

For better visualization of vibrational modes, 
the mode shapes are superimposed on a representa- 
tion of sun, planet and ring gear base circles. 

The equilibrium position Is represented by the line 
of action of each gear mesh between points of tan- 
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gency on respective base circles^ and by radii 
through these points of tangency. Rotational dis* 
placements of these tangency points (computed as 
Ra displacements along the line of action) define 
angular displacements of the planet. Any transla- 
tion Is added to the computed rotational displace- 
ments. The equilbrium position of a single planet 
at arbitrary angle a Is shown in Table 2(a) and 
Illustrated In Fig. 7. Deflections from equilibrium 
due to rotation are shown In Table 2(b). 

Verification of Nethod 

To validate the algorithm presented above, 
results were compared with those of earlier studies. 
Cunllffe^ performed analytical and experimental 
studies on a small star gear box with the ring gear 
rigidly fixed In translation and the carrier fixed 
In both rotation end translation, thus eliminating 
5 degrees of freedom from the 18 degrees of freedom 
model. Using the Cunllffe data (suppressing the 
unwanted modes with extremely high - - stiff- 

ness values), 5 axisymmetric modes and 4 orthogonal 
pairs of asyiMtrIc modes were duplicated* As an 
example, the highest (torsional) mode (which ex- 
hibits anti synchronous motion at the sun/planet 
mesh) Is shown In Fig* 8. 

Botman^ Investigated an 18 degrees of free- 
dom planetary gear system to obtain 6 symmetric 
modes and 6 pairs of asymmetric modes. Similar 
results are obtained In the present work. The 18 
modes are described in Table 3. Two of the highest 
modes obtained with the Botman model are presented 
In Figs. 9(a) and 10(a). These are very similar to 
the Cunllffe 13 degrees of freedom torsional modes 
despite differences In relative sizes of the compo- 
nents and In the additional degrees of freedom. 

Results and Discussion 

The algorithm developed In this Investigation 
was used to perform a series of comparative parame- 
tric studies to assess the effect of unequal mesh 
stiffnesses on the natural frequencies and the re- 
spective mode shapes for a planetary gear system. 
Three cases are reported: (1) All six mesh stiff- 
nesses equal; (2) One sun/planet stiffness reduced 
(to one half the nominal value); (3) One planet/ring 
mesh stiffness reduced. Brief consideration was 
also given to simultaneous stiffness reductions at 
two meshes to simulate variation In nun4>er of teeth 
In contact. 

Figure 9 shows the highest axlsymmet'^lc mode 
(ho. 18 In Table 3) which exhibits ant1-:,ynchronous 
motion of the planets with respect to the sun and 
ring gears. In Fig. 9(a) all mesh stiffnesses are 
equal (case 1). Figure 9(b) Illustrates the effect 
of reducing one sun/planet mesh stiffness (case 2). 
Figure 9(c) shows the effect of reducing one 
planet/ring mesh stiffness (case 3). Figure 10 
Illustrates cases 1, 2, and 3 respectively for one 
of the orthogonal pair of related asymmetric modes 
(No. 16 In Table 3). Figure 10(b) Is particularly 
affected by the reduced stiffness at the sun mesh 
which allows significant lateral motion of one 
planet. 

Although some mode shapes changed noticeably 
due to unequal stiffnesses, the values of the 
natural frequencies were little affected, varying 
less than 5 percent for a 50 percent reduction In 
the stiffness at one mesh. This Indicates that the 
eigenvalues computed from average stiffnesses are a 
good approximation for the forced response solution 


of a 4ytimic system which has time varyinq stiff* 
nesses due to the alternation between one and two 
pairs of teeth In mesh. 

Natural frequencies were used to determine the 
time step and total Integration time for the numer- 
ical calculation of dynamic load factors (Eq. (1) 
with time-varying, periodic stiffnesses and without 
the assumption of zero damping). The dynamic load 
factor Is defined as the ratio between the maximum 
Instantaneous load and the static transmitted 
(design) load. Five percent of the system's short- 
est natural period provided a stable time step for 
the integration process. To assure the decay of 
start-up transients, the Integration was carried 
out for ten times the longest natural period. 

Figure 11 shows the variation In the dynamic 
load factor for one sun/planet and one planet/ring 
mesh engagement for the planetary gear system out- 
lined In Table 4. The peak dynamic load factor 
occurs near the frequency of the hl^st mode found 
by the eigenvalue method outlined above and Is 
excited by the tooth meshing frequency. Identifi- 
cation of this resonance frequency Is particularly 
Important since (for the design analyzed here) 
dynmnlc factors above 2.5 are Indicative of tooth 
separation, and constitute an unacceptable 
operating condition. 

Summary of Results 

An algorithm suitable for a mini-computer was 
developed for finding the natural frequencies and 
mode shapes of a planetary gear system which has 
unequal stiffnesses between the sun/planet and 
planet/ring gear meshes. Mode shapes are re- 
presented In the form of graphical computer output 
that illustrates the lateral and rotational motion 
of the three co-axial gears and the planet gears* 
This procedure permits the analysis of gear trains 
utilizing non-uniform mesh conditions and user 
specified masses, stiffnesses, and boundary condi- 
tions. Numerical Integration of the equations of 
motion for a three planetary gear systems was per- 
formed. The following results were obtained: 

1. Excellent agreement was obtained between 
natural frequencies calculated assuming constant 
mesh stiffnesses and resonant frequencies of the 
dynamic system. 

2. The algorithm can be used for determining 
the time step and total Integration times In the 
numerical solution of the differential equations of 
motion. 

3. The methodology Is a quick and efficient 
means of Identifying operating resonance regions 
and analytically tuning the drive train by examining 
the effect of parameters such as planet and sun gear 
support stiffnesses. Input /output shaft stiffnesses, 
and gear mesh stiffnesses. 

4. The algorithm Is adaptable for use In a 
Computer-Aided Design (CAO) environment. 
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TABLE 1. - STIFPNESS TERNS 
(a) Bearing support reactions. 


Component Displacement T/q F^/q Fy/q 


Sun 

Planet-J 

Carrier 

Ring 



(b) Sun/planet Interaction. 



T/q 

Fx/fl 

-KSP Rj 

-KSP COS a 

-KSP Rj cos ft 

-KSP COS^« 

-KSP sin ft 

-KSP cos ft sin ft 

-KSP Rp 

-KSP COS a 

-KSP Rp cos ft 

-KSP COS^B 

-KSP K « 

-KSP cos 0 sin ft 

p 



Fy/q 


-KSP sin a 
-KSP cos ft sin ft 


-KSP sin ft 
-KSP cos ft sin ft 


(c) Carrier/planet Interaction. 


Component I Displacement I T/q 



Planet "j" Xj - X 


(d) Ring/planet interaction. 
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TABLE 2. - NODE SHAKS 

(a) Egulllbrluai position: center and base circle tangency points. 




TABLE 3. - DESCRIPTION OF MODE SHAPES 


Mode no. 

Type 

Motion 

Strain energy 

lb 

3 

Translational 

Symmetric 

Sun 

Sun mesh 

4 

5a, 5b 

Symmetric 

Translational 

Ring, carrier 

Ring mesh, output shaft 

7a, 7b 
9 

Translational 

Symmetric 

Radial planet motion 

Planet bearings 

10 

11a, 11b 

Symmetric 

Translational 

Carrier, sun 

Output shaft 

13a, 13b 
15 

Translational 

Symmetric 

Tangential planet motion 

Ring mesh, planet bearings 

16a. 16b 
18 

Translational i 
Symmetric 

Anti synchronous motion 

Sun mesh, ring mesh 
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TABLE 4. - DESIGN AND SYSTEM PARAMETERS USED 
FDR CALCULATING DYNAMIC FACTORS 


Design Parameters 

Diametral pitch 5 

Pressure angle, deg 22,5 

Number of teeth 

Sun 25 

Planet 24 

Ring 73 

Face width. In 1 

Theoretical contact ratio 

Sun/planet 1*56 

Planet/ring .... 1*78 

System Properties 

Jdriver. in-lb«sec2 5.M0 

Jjun* in-lb* sec^ 0.040 

Msun* lb*sec2/1n 0.015 

JolSnet. in-Ib-sec? 0.033 

Jhw. 2.000 

Jload. 1n-lb»sec2 5.000 
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(a) Partial cutaway. 

(b) Detail, planet carrier assembly. 
Figure 1. - Planetary gear train. 
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Figure 2. - Planetary gear stage showing global and local coordinate 
systems The direction of positive rotation for each component is 
indicated by arrows. 
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Figure 3. - Dynamic modei of pianatary grer system. 
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Figure 5. - Carrier-planet interaction for planet "j" at arbitrary 

angles 


RMG 

GEAR 


^fGfNAL PAG- 
^ POOR quality 




Figure 7. - Equilit)rium position of pianot "J" at arbitrary angla Tha 
direction of positive rotation for sun. pianet and ring ^ars is indicated 
by arrows. 
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Figure 1 - 13 OOF Antisynchronous $ym* 
metric mode. 
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(a) All mesh stiffnesses equal. 

(b) One sun/pianet mesh stiffness reduced. 
Ic) One pianet/ring n.esh stiffness reduced. 

Figure 9. - 18 DOF Ai.dsynchronous sym- 
metric mode. 


(a) All mesh stiffnesses equal. 

(b) One sun/planet mesh stiffness reduced. 

(c) One planet/ring mesh stiffness reduced. 

Figure 10. - 18 DOF Antisynchronous tran- 
slational mode. 
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0 1 2 3 4 3 6 7 8 9 10 11 12 13 14 15 16 

INPUT SPEED, 1000 rpm 


(al Sun/planet mesh. 

(b» Planet/ring mesh . 

Figure 11. - Dynamic factor as a function of sun input speed. 


